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Abstract. For the electromagnetic interaction of two particles the relativistic quantum 
mechanics equations are proposed. These equations are solved for the case when one particle 
has a small mass and moves freely. The initial wave functions are supposed to be concentrated 
at the coordinates origin. The energy spectrum of another particle wave function is defined 
by the initial wave function of the free moving particle. Choosing the initial wave functions 
of the free moving particle it is possible to obtain a practically arbitrary energy spectrum. 



fi; 1 Introduction 

For the particles of the mass \x > and the spin 1/2 the relativistic equation was proposed 
ON by Dirac ([1], equation (4.29); [2], equation (1 - 41); [3], equation (2.9)) 

1 i n 



4 3 / 

E E (7*)nm - 
n=lu=0 \ 



, i— — \xl) n (x)+ ni) m (x) = Q. (1.1) 



Let us consider a classical electromagnetic field with a vector potential A v (x). If we change in 
the equation (1.1) the differential operator —id v for the differential operator —id„ + qA v (x), 
we get the equation for the interaction between the particle with the charge q and the 
external classical electromagnetic field ([1], Chapter 4, equation (4.202); [3], Chapter 2, 
equation (2.62)) 

43 / d \ 

E E (VOmn -i^TT, + 1 A »( X ) 4>n(x) + ^m{x) = 0. (1.2) 

The substitution of the electron charge q = — e and the Coulomb vector potential A (x) = 
Ze(47r|x|) _1 , A k (x) = 0, k = 1,2,3, into the equation (1.2) yields the equation for the 
electron in the external electromagnetic field generated by the nucleus with the charge Ze. 
We consider a nucleus as a classical particle. The energy spectrum of this equation (1.2) is 
discrete and the energy level E 1>1/2 = /uc 2 (l - (Z/137) 2 ) 1/2 ([3], equation (2.87)). For Z > 
137 the value £7i,i/2 becomes imaginary. The description of the electromagnetic interaction 
between the electron and the nucleus with the charge Ze > 137e seems not to make sense. 
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At present time the nuclei with the charges Ze < 118e are synthesized. (The nucleus 
with the charge Ze = 117e has not been yet discovered.) The equation (1.2) should be 
mathematically self - consistent for Coulomb vector potential with any charge. Due to 
the paper [4] Hamiltonian (1.2) with the Coulomb vector potential A (x) = Ze(47r|x|) _1 , 
A k (x) = 0, k = 1,2,3, has a self adjoint operator extension for any charge Ze. For Ze > 
(y / 3/2)137e > 118e this self adjoint operator extension is not unique. 

The equation (1.2) defines the wave function of the electron interacting with the elec- 
tromagnetic field generated by the classical particle. Both interacting particles should be 
quantum. In the quantum electrodynamics ([5], Lecture 24) the following equations for the 
electromagnetic interaction of two particles are studied 

E f fll fS(7")m.,n. [ijrz] -VsSm s ,nX) ^ ni n 2 , PlP2 (xi, X 2 ) + ]T 
ni ,n 2 = l\\s = l V = \ UX sJ J J ui,v 2 = 

■q^Kq^DKxi - x 2 ) J| (7^) m s ,n s I ' l Pnin2,p\P2 

(Xi, X 2 )\ =«[[ $( X s)8m s ,p s , (1.3) 

D c m2 (x) = £ lim o (27r)- 4 J d 4 kexp{i(k,x)}(m 2 - (Jfe, Jfe) - ie) _1 , 

(x,y) = x°y°- X>V, (1-4) 

k=l 

where q 2 , q 2 are the charges, K is the electromagnetic interaction constant, r/^ = rf v is the 
diagonal matrix with diagonal matrix elements i]q = —rju = —r\ 22 = —i] 33 = 1. 

The solutions of the equations (1.3) contain the divergent integrals [5]. The finite answers 
for the divergent integrals are obtained by means of the renormalization procedure. Due to 
the book ([6], Chapter 4): "The shell game that we play to find n and j is technically called 
"renormalization". But no matter how clever the word, it is what I would call a dippy 
process! Having to resort to such hocus - pocus has prevented us from proving that the 
theory of quantum electrodynamics is mathematically self - consistent. It's surprising that 
the theory still hasn't been proved self - consistent one way or the other by now; I suspect 
that renormalization is not mathematically legitimate. What is certain is that we do not 
have a good mathematical way to describe the theory of quantum electrodynamics: such 
a bunch of words to describe the connection between n and j and m and e is not good 
mathematics." 

The equations (1.3) do not satisfy the causality condition: the support of the fundamental 
solution Dq(x) of the wave equation does not lie in the upper or lower light cones. The distri- 
bution Dq(x) choice is probably connected with the causality condition. The chronological 
product of two free scalar field operators is defined as 

T(d>(x)(f>(y)) = (f>{x)(t>{y) + ie m , (y - x)I (1.5) 

where I is the identity operator and the distribution 

-e m 2(-x) = lim^Y 4 j d 4 kexp{i(k,x)}(m 2 - {k° + ief + 1 - (1.6) 

The vacuum expectation of the product of two free scalar fields is defined as 

< (/>(x)(/>(y) > = -iD m2 (x - y), 
D~ 2 {x) = i{2n)- 3 J rf 4 A;exp{-i(A;,a:)}^(A; )5((A;,A;) -m 2 ). (1.7) 



2 



The vacuum expectation of the chronological product (1.5) is equal to 



< T{4>{x)4>{y)) > = -iD~ 2 (x - y) + ie m ,{y - x) = -iD c ^{x - y). (1.8) 

Stueckelberg and Rivier [7] believed that the classical " causal action" is given by the distri- 
bution — e m 2 {y — x) and the distribution (x — y) = D c m 2 (y — x) defines the probability 
amplitude of the "causal action". 

The equations (1.3) do not look like the equations (1.2). We want to write down two 
equations of the type (1.2) where one particle interacts with the electromagnetic field gen- 
erated by another particle. These equations should satisfy the causality condition. Let us 
define 4x4- matrices 

■w-('riH-(^)^-(^ , ')."-« «-> 

where 2x2- matrices a u are given by the relations 

•■-(; ;)■■•-(! ;)■••-(! o)^-(;-i)- 

Let us consider two equations for the wave functions (ip s )n s ( x s), s — 1)2, equal to zero for 
the negative x° s 

4 

-i E (7°)m.,n.5(a:2)(^)n.(+0,x i ), x s = (xl,x 2 s ,x s s ) E R 3 , s = 1,2, (1.11) 

n 3 = 1 

where (ip s ) ns (+0, x s ) are the initial wave functions. For //i,/i2 > the equations (1.11) with 
the right - hand sides equal to zero are equivalent to Dirac equation (1.1). For /ii,//2 > 
the equations (1.11) with the right - hand sides equal to zero are equivalent to the system 
of Klein - Gordon equations. The right - hand sides of the equations (1.11) are the standard 
method to define Cauchy problem. The equations (1.11) describe the free motion of particles 
with the masses 1^1,1^2 > 0. In contrast with the equation (1.1) the particles with zero mass 
are not distinguished. If the initial wave functions (0s)n s (+O, x s ) = (ip s ) ns 5(x s ), s = 1,2, 
then the product of the equations (1.11) looks likes the equation (1.3) with the interaction 
constant K — 0. 

The equation of the electromagnetic action of the second particle on the first particle has 
the form 

E 4 / ^(il {ilyUn, + (/«))~) Mn s (Xs) + 



3 = 1,2 



3 2 n 

Kq^Tj^eoix! - X 2 ) J[ (Y S )m 3 n s Mn s ( X s) 
ui,V2=0 s = l j 

2 



E / II (tV^X^W+O.X.). (1.12) 

= l,-,4, J s = l 



n s 

s = 1,2 
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The equation of the electromagnetic action of the first particle on the second particle is the 
similar one. Both equations of the type (1.2) satisfy the causality condition. The support 
of the distribution — e ( — x) lies in the closed lower light cone. The fundamental solution 
— e (— x) of the wave equation is unique in the class of the distributions with supports in the 
lower light cone. The causality condition defines the distribution — e (— x) uniquely. The 
distribution (1.6) differs from the distribution (1.4) in the rule of going around the poles in 
the integral. The support of the distribution — eo(— x) lies on the lower light cone boundary. 
The distribution — e (— x) gives the delay. It is necessary to have the delay in a relativistic 
interaction equation. It was already noted by Poincare [8]. The interaction propagates not 
instantly but at the speed of light. We have to take into account the distance between the 
interacting particles. The delay is one of possible causality statements. 

In this paper the equations (1.12) are solved for the case when the second particle has 
a small mass and the wave function (ip 2 ) n2 {.x 2 ) satisfies the second equation (1.11) for the 
initial wave function (ip 2 ) n2 (+0, x 2) = (V , 2)n 2 ^( x 2)- The energy spectrum of the solutions 
ijl 'i)m( x i) is defined by the initial wave function {i^ 2 ) n2 5{yi 2 ). By making a choice of the 
vectors (ip 2 ) n2 it is possible to obtain a practically arbitrary energy spectrum. 



2 Dirac equation 

It is necessary to write down an equation of the type (1.1) for the free motion of a particle 
with arbitrary spin and mass. Let us consider the complex 2x2- matrices 



The 2 x 2 - matrix 

/ A.. A„ \ 

(2.2) 



A u A 21 

A\ 2 A 22 



is called Hermitian conjugate. If A* = A, the matrix (2.1) is Hermitian. The matrices (1.10) 
form a basis of Hermitian matrices. The multiplication rules for the matrices (1.10) are 

trV = a , a°a^ = a»a° = <t", fi = 0, 3; 

3 

a k 1(T k 2 = J- e klk2k Ha ks , h,k 2 = 1,2,3, h ^ k 2 (2.3) 
fe 3 = i 

where the antisymmetric tensor e klk2ks has the normalization e 123 = I. 

We identify the four dimensional Minkowski space with the four dimensional space of 
Hermitian 2x2- matrices 

3 

x = xW- (2-4) 

/i=0 

For a complex 2x2- matrix (2.1) we define the following 2x2- matrices 

a t =( A a u ! 2I V > f ! ;i t 2 V (2-5) 

V A 12 A 22 J ' ^ A 21 A 22 j K ' 

The matrices (2.1) with determinant equal to 1 form the group SX(2, C). The matrices 
(2.1) satisfying the equations A* A = a , det A = 1 form the group SU(2). The group 
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577(2) is the maximal compact subgroup of the group SL(2, C). Let us describe the irre- 
ducible representations of the group SU(2). We consider the non - negative half - integers 
/ = 0, 1/2, 1,3/2, .... We define the representation of the group SU{2) on the space of the 
polynomials with degrees less than or equal to 21 

TtAMz) = (A 12 z + /4 22 f ( ^ + ^ ) . (2.6) 

We consider a half - integer n = —/,—/ + 1, — 1, 1 and choose the polynomial basis 

Mz) = ((l-n)\(l + n)\)- 1 / 2 z l - n . (2.7) 
The definitions (2.6), (2.7) imply 

T,(Aty B (z) = £ i; m (z)t l mn (A), (2.8) 

m=— I 



tLW = ((I ~ m)\{l + m)\{l - n)\{l + n)\)^ x 

oo A-m-j aj A m-n+j *l+n-j 

\ - ^11 /i 12 /1 21 ^22 /<-, qx 

■f^ r(j + l)r(Z - m - j + l)T{m-n + j + l)T(l + n- j + 1 ' j 

where T(z) is the gamma - function. The function (r(^)) -1 equals zero for z — 0, —1, —2, .... 
Therefore the series (2.9) is the polynomial. 

The relation (2.6) defines the representation of the group SU(2). Thus the polynomial 
(2.9) defines the representation of the group SU(2) 

t l mn (AB) = £ t l mk (A)t l kn (B). (2.10) 

k=-l 

This {21 + 1) - dimensional representation is irreducible ([9], Chapter III, Section 2.3). The 
relations (2.9), (2.10) have an analytic continuation to all matrices (2.1). 

By making the change j^j + n — moi the summation variable in the equality (2.9) we 
have 

t l mn {A)=t l nm {A T ). (2.11) 

The definition (2.9) implies 

= U (2.12) 

The polynomial (2.9) is homogeneous of the matrix elements (2.1). Its degree is 21. The 
sum (2.3) contains the only non - zero term. The relations (2.3), (2.10) imply 

E t E.(ECK)4 P (^) (-^)) (£/^n(Oi(^) (-^)) = 

E 4n(^ fc3 )in(-^)(( e<!lfc2fc3 ) 2 ' +2/ + (e k2klki ) 2l+21 )- ~ 



l<fc 1 <fc 2 <3, UJj udj 

fc 3 = l,...,3 



'^mn^mh{d x , d x ) , (2.13) 
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(8.,8.)=T,V m (J^ ■ (2-14) 
For an odd integer 21 + 21 the relation (2.13) has the form 

E E fE4 P (04(^)f-^)VEoU04(^)^- 9 v ' 



p = -lp=-l \v = \ U It \v = \ ux II 

-8 mn 5rhh(d x ,d x ). (2.15) 

By making use of the relation (2.15) for an odd integer 21 + 21 Lorentz invariant equations 

(-(d x ,d x )-fi 2 ) ( j )m ^(x)=0, (2.16) 

m = —I, —I + 1, I — 1, Z, m — —I, + — may be rewritten as the system of the 

linear equations 

4Z+2 2/+1 3 / o \ 

E E E »? w 4-i-i,n-3i-2(04-i-i lA -i-i(^) <M*) + 

n = 21+2 h = lv = \ ux I 

ip m m{x) = 0, m = 1, 2/ + 1, m = 1, 21 + 1; 

2Z+1 2/+1 3 / Q \ 

E E E 4-3i-2,„-,-i(04-i-i,*-i-i(^) iM*) + 

n = ln=li/ = \ ux I 

/iV m mW=0, m = 2/ + 2,...,4/ + 2, m = 1, 21 + 1. (2.17) 
Let us define the ((4/ + 2)(2Z + 1)) x ((4/ + 2)(2Z + 1)) - matrices 

(<Xl,i(t* 2 ))mm,nn = l^^rnJfnh, m, 77 = 1, 21 + 1, 777, 77 = 1, 2Z + 1; 

(a^i/J 2 ))^^ = 5 mn 5rnn, 771, 77 = 2/ + 2, . . . , 4/ + 2, 777, 77 = 1, 2Z + 1; 

(Pl,i(V 2 ))mm,nn = $mn5rhh, 777, 77 = 1, 2/ + 1, 777, 77 = 1 , . . . , 2Z + 1; 

(Plj(^ 2 ))mm,nn = ^^mn^rhh, 777, 77 = 2/ + 2, . . . , 4/ + 2, 777, 77 = 1, . . . , 2Z + 1; 

(lli( a °))mrh,nn = r f V ^m-l~\,n-Zl-2 {^^m-i-lfi-l-l ' 
777 = 1,...,2Z+ 1, 77 = 2/ + 2, ...,4/ + 2, 777,77 = 1,..., 2/ + 1, 1/ = 0, 3; 



(7n - ( <7 °))"»»">i" ~~ tm-3l-2,n-l-l( a K 



m— Z— l,n— /— 1 



a" 



777 = 21 + 2,..., 4/ + 2, 77 = 1,...,2/ + 1, 777,77 = 1,..., 2/+ 1, 1/ = 0, 3. (2.18) 

The other matrix elements are equal to zero. If an integer 2Z + 21 is odd, then an integer 
(2Z + 1)(2Z + 1) = 4ZZ + 2Z + 2Z + 1 is even and 4 divides into (4Z + 2)(2Z + 1). For 1 = 1/2,1 = 
the integer (4Z + 2)(2Z + 1) = 4. The definitions (2.18) are the straightforward generalizations 
of the definitions (1.9). The definition (2.9) implies 

tL(A) = A m+ s tn+ s. (2.19) 

4 x 4 - matrices ai/2,o(^ 2 ), Pi/2,o(^ 2 ), 7i/2,o( cr °) coincide with the matrices a(/7 2 ), /3(/7 2 ), 
Y given by the relations (1.9). By making use of the definitions (2.18) we can rewrite the 
equations (2.17) as 

Al+2 21+1 / 3 ( d \ \ 

E E E (Tf/O Wi ^T— + (A,K^))mrh,nn = 0. (2.20) 



n=lti=l \i/ = 
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In contrast with the equation (1.1) the mass \i in the equation (2.20) may be equal to zero. 
The definitions (2.18) imply 

M^)A> 2 )=M,M (2-21) 

In view of the relations (2.21) the action of the matrix a l [(/j.) on the equation (2.20) yields 
the equation of type (1.1) 

41+2 21+1 3 / g \ 

E E E Wj(<r°))mm, n n ( J ^n(x) + fiUrn(x) = 0, (2.22) 



n = l h = l v = 



41+2 21+1 

irnfn{x) = E E (Plj(f^)) mrh,nh'4'nh(x) = 0. (2.23) 
n = 1 h = 1 

For n > the transformation given by the relation (2.23) is the isomorphism. Due to the 
relations (1.9), (2.18), (2.19) the equation (2.22) for I = 1/2, Z = coincides with Dirac 
equation ([2], equation (1 - 41)). The relations (2.18) imply 

{a lj i{n 2 )P lt i{n 2 ))mm,nh = {Pl^ 2 )^l^ 2 ))mrh,nh = l??>mJmhi (2.24) 

to, n — 1, 4/ + 2, to, n = 1, 21 + I. In view of the second relation (2.21) and the relation 
(2.24) the action of the matrix 7^(<t°)q; ; [(fi 2 ) on the equation (2.20) yields 

41+2 21+1 / 3 / Q \ \ 

E E (E(7?^ )7M(^)A,K^)) mm ,nn \-iQ^J +^(liy)) m m,nnj ^nn(x) = 0. 

(2.25) 

The relations (2.3), (2.10), (2.12), (2.18) imply 

((7°/(° r °)) 2 )mm,nr i = ^mn^mn, 771,71=1, 4/ + 2, TO, 77, = 1, 21 + 1, (2.26) 

(7/ ; /( cr )T; i ;'( cr ))mih,nh (^i /(^ )7/ ) /'( cr ))mm,nn 

tm-Z-l,n-i-l (°" fc ) t !„_/_i in _/_i ) 

771, n = 1, 2/ + 1, m, 72 = 1, 2/ + 1, A: = 1, 2, 3; 

(^,/( <7 )^,/( <J ))i™h,nh — ~(ll j i( <J )7/ ) /(°" ))mm,nh = 
- t m-3l-2,n-3l-2( Crk ) tl m -i-l,n-i~l(°' k )i 

m, n = 21 + 2, 4/ + 2, to, n = 1, 2Z + 1, fc = 1, 2, 3. (2.27) 

The other matrix elements are equal to zero. The coefficients of the polynomial (2.9) are 
real. Hence in view of the relations (2.11), (2.18), (2.26), (2.27) the matrices 7,°,-(<7°), 
/ y^(a°)^:(a°)f3 l i(n 2 ), v — 0, ...,3, are Hermitian. Due to the relation (2.26) we have 
(7^(<7°)) 2 /3 M -(/z 2 ) = /3 M -(/i 2 ). Now the equation (2.25) implies that the integral 



41+2 21+1 

/d 3 x£ E(A,K/^ 2 ))™— \i>mrn{x)\ 2 (2.28) 
m = l ih = l 
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is independent of the variable x° for x° > 0. The integrand (2.28) is called the probability 
density of a solution of the equation (2.20). For ii > the transformation given by the 
relation (2.23) is the isomorphism. The probability density (2.28) expressed through the 
function (2.23) is independent of the mass and coincides with the usual probability density 
for the function (2.23). In the quantum mechanics the fixed probability density defines 
Hilbert space where any Hamiltonian acts. The integral (2.28) depends on the parameter \j? 
in the equation (2.20). We do not want to suppose that all particles have strictly positive 
masses. The experiments deal with the asymptotic solutions of the interaction equations. We 
expect the solutions of the interaction equations coincide asymptotically with the products 
of the equation (2.20) solutions. The probability density of the last solutions is given by the 
integrand (2.28). 

The relations (2.3), (2.10), (2.12), (2.18), (2.26), (2.27) imply 

(^}(a )^(a )) mrh>nfl = -5 klk2 6 mn S mn + £ {-if x 

k 3 = i 

((a Zi K0)7 / ^ )73(^°)) mm ,n. - (^/(OJtJ/MtJM)™*,^) > 

m, n = 1, 41 + 4, rh, h = 1, 21 + 1, k 1: k 2 = 1, 2, 3. (2.29) 

For an odd integer 21 + 21 the relations (2.26), (2.27), (2.29) imply 

\ lii\ u ) ) )mm,nn * \ lljy ) ";/V ) )mm,nn ^'1 u mm,nni 

pi, v = 0, 3, m, n = 1, 4Z + 2, rh, h = 1, 21 + 1. (2.30) 
Let us define the representation of the group SL(2, C) in the Lorentz group 

3 

]T A$(A)x v <T' i = AxA*. (2.31) 

We define also the (4/ + 2) (2/ + 1) - dimensional representation of the group SL(2, C) 

(A), 

m, n — 1, 21 + 1, rh, h — 1, 21 + 1; 

rh—l—l,n—l—l 

m, n = 21 + 2, 41 + 2, rh,h=l, 21 + 1. (2.32) 

The other matrix elements are equal to zero. The definitions (2.18), (2.32) imply 

StMfaMWM' 1 ) = MA ^(A^ji^S^A- 1 ) = /3 M (/i 2 ). (2.33) 

Let the functions ip m m{x), m — 1, ■■■,41 + 2, rh = 1, 21 + 1 be the solutions of the equation 
(2.20). The relations (2.33) imply that the functions 

41+2 2/+1 

Um{x) = E E (S^AyU^nniA-^iAT 1 ) (2.34) 

n= 1 h= 1 
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are the solutions of the equation 



41+2 21+1 / 3 / a \ \ 

E E (E (7j(^)U,»« ("WJ + (^^ 2 ))™*.»*J W*) = °' ( 2 - 35 ) 
7^) = E KimMK^iM' 1 ) ( 2 - 36 ) 



i/ = 

for any matrix A G SL(2, C). The definition (2.36) implies 

3 

z 

u = 



= E A^(A)5 M -(A) 7 ^.(5)5 M -(A- 1 ) (2.37) 



for any matrices A,B G SL(2, C). By changing the coordinate system we change the matrix 
7^-(<7°) in the equation (2.20) for the matrix (2.36). The solutions of the equation (2.20) are 
transformed to the solutions (2.34) of the equation (2.35). It is valid for all non - negative 
half - integers I, I. Due to the book ([2], relation (1 - 43)) for any matrix A G SL(2, C) 



7^)=^ (O (2-38) 

2 1 



Hence the equation (2.20) for I — 1/2, I — is covariant under the group SL(2, C). The 
relation (2.38) picks out Dirac equation. 

For an odd integer 21 + 21 the relations (2.30), (2.36) imply 

fi, v = 0, 3, to, n = 1, 41 + 2, to, n = 1, 2Z + 1 (2.39) 

for any matrix A G SX(2, C). 

Let the functions (, m rh(x) be the solutions of the equation (2.35). Let us introduce the 
distributions 

fmrn(x) = 9( X °)U m (x), 0(x) = j J' ^ jj' (2.40) 

The equation (2.35) implies 

41+2 21+1 / 3 f 8 \ \ 

E E ( E (7n'(^))mm,Bn \ ) )) m " 1 .™™) fnn{ x ) = 

-«(x°)/!L(+0,x), (2.41) 

4i+2 2Z+1 

/^(+0,X)= £ E(7°/(^))mm,nn/nn(+0,x). 
n = 1 n = 1 

Let a distribution ^2 (re) G S"(R 4 ) with support in the closed upper light cone satisfy 
the equation 

(iih^a) -^)) ^i,..A x ) = *(*)■ ( 2 - 42 ) 

We prove the uniqueness of the equation (2.42) solution in the class of the distributions 
with supports in the closed upper light cone. Let the equation (2.42) have two solutions 
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e i 2(^)1 e 2 2 (^)- Since its supports lie in the closed upper light cone the convolution 



is defined. Now the convolution commutativity implies these distribution coincidence: 



n \ r 

n H4A) -/*?)/ <*vsu<* - = «suw- (2-43) 



Vi = l 



Due to the book ([10], Section 30) 

e (x) = - (2n)- 1 9(x )5((x, x)), e 0fi (x) = (87r)-^(x )^((x, a;)). (2.44) 

The uniqueness of the equation (2.42) solution in the class of the distributions with supports 
in the closed upper light cone implies 

-(d x ,d x )eo fi (x) = e (x). (2.45) 

In view of the second definition (2.40) 

(d x ,d x )(9(x°)9((x,x))(x,x) n ) = 4n(n + l)9(x°)9((x,x))(x,x) n - 1 , n = 1,2,.... (2.46) 

Due to the relations (2.44), (2.46) the distribution e ,...,o(^) with n zeros has the form 

e ,...,o(z) = (-l) n (27r4^ 1 (n - 2)!(n - l)!)- 1 ^ )^, a;)"" 2 , n = 2, 3, .... (2.47) 

Let us prove the following equality 

e^(x) = - (27r)- 1 9(x°)5((x,x)) + 

00 

9(x°)9((x, x)) ]T /if t (-l) n+1 (27r4 n (n - l)!™!)" 1 ^, x) n -\ (2.48) 

n = l 

By making use of the relations (2.44), (2.45), (2.46) it is possible to prove that the right - 
hand side of the equality (2.48) satisfies the equation (2.42) for n = 1. The support of the 
right - hand side of the equality (2.48) lies in the closed upper light cone. The solution of 
the equation (2.42) is unique in the class of the distributions with supports in the closed 
upper light cone. 
Let us prove 



d 4 pexp{-i(p,x)} Hdp' + ie) 2 -^-^)- 1 . (2.49) 
j = i 

The integral (2.49) is the solution of the equation (2.42). By making the shift of the integra- 
tion path in the right - hand side of the equality (2.49) we obtain that the distribution (2.49) 
is equal to zero for x° < 0. The distribution (2.49) is Lorentz invariant. Hence its support 
lies in the closed upper light cone. Now the uniqueness of the distribution e^j ^2 (x) implies 
the equality (2.49). The equality (1.6) is the particular case of the equality (2.49). 
The relations (2.21), (2.24), (2.39) imply the relation 



4/+2 2/+1 / 3 / a \ 

EE E (7^))™ H^; + 

p=\p=\ \ v = \ / 



mm,pp I X 



\u = 



dx u , 

m, n = 1, 4Z + 2, m, h = 1, 2Z + 1, (2.50) 
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for an odd integer 21 + 21. The relations (2.42), (2.50) imply that the solution of the equation 
(2.41) has the form 

4Z+2 2/+1 / 3 f 8 \ \ 

fmrh{ x ) — E E ( E (7i ; /(^))mrh,nn I ~ ^~o37 J ~~ { a l,ii^ ))mm,nh ] x 
ji = ln = l\i/ = \ OX J J 

(-i J d 4 ye^(x-y)5(y )f^(+0,y^ , (2.51) 

m = 1, ...,4/ + 2, m — 1, ...,2/ + 1. For the solution of the equation (2.35) in the domain 
x° < it is sufficient to use the distribution — e^(— x) in the relation (2.51). 

We suppose that the smooth function / n ^(+0, x) is rapidly decreasing at the infinity. By 
shifting the integration path in the integral (2.49) we have 

J d i ye lt ,(x-y)6(y )f n *(+0,y) = 

(27t)- 4 | d* P eMx -KP,x)}((P° + *) 2 - |p| 2 -^r7nn(+0,-)(p), 

r 3 

/mm(+0,-)(p) = / rf 3 xexp{-2^p A; x fc }/ mril (+0,x). (2.52) 

J k=l 

The integral with respect to p° may be easily calculated. For x° > and / mr h(+0, -)(p) = 
/mm^p-q) the functions (2.51), (2.52) are not the eigenfunctions of the differential operator 
—id/dx° and are the eigenfunctions of the differential operator (—id/dx ) 2 . 



3 Relativistic quantum Coulomb law 

We consider at first the relativistic Coulomb law in the classical mechanics. The relativistic 
Coulomb law is the particular case of the relativistic Newton second law 



ds dt \ ds dt 



ds dt 



dtd dtdx^\ q Jt, UUT , , ,dtdx ai 

>nc^^ I - — - I + - E E vTF IMi ... ak {x)- 

C k = 0ai,...,a k =0 

dt 
ds 



dt dx ak 
ds dt 



1 ( 1 - c- 2 






dt 



= o, 



-1/2 



(3.1) 



where x° — ct, /j, — 0, ...,3. In the equation (3.1) the force is the polynomial of the speed. 
For an infinite series of the speed it is necessary to define the series convergence. The second 
relation (3.1) implies the identities 



Vota 

a = 



dt dx a ^ 



3 dtdx a dt d (dtdx aS 

1, 2 , '/<■" 

a = 



ds dt I r n ds dt ds dt \ ds dt 

/ a = u \ 

The equation (3.1) and the second identity (3.2) imply 

, , dt dx ai dt dx ak +^ 



0. 



(3.2) 



N 



E E 

k = ai,...,a k+1 =0 



ds dt 



0. 



(3.3) 
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Let the functions F ai ... a (x) satisfy the equation (3.3). Then three equations (3.1) for 
\i = 1, 2, 3 are independent 

j N 

* ((1 _ c -2| v | 2) -l/2^ _ g £ c -*(l - C - 2 |v| 2 )- V X 

fe = 

F te ...„ >W =0, «• = _,, = 1,2, 3. (3,4) 

<*i,...,a fc = U 

The following lemma is proved in the paper [11]. 

Lemma. Let there exist a Lagrange function L(x, v,t) such that for any world line, x^it), 
x°(t) = ct, and for any i = 1,2,3 the relation 

4™-™=* ((1 - c- 2 |v| 2 )-VV) - 
dt dv l dx l dt V v 



x 3 dx ai d.r" 1 



gE^(l-c- 2 |v| 2 )-^ E ~TT~ (3-5) 

fc = ai,...,Qfc = 

holds. Then the Lagrange function has the form 

n 3 

L(x, v, t) = -mc 2 (l - c~ 2 |v| 2 ) 1/2 + -Y, ty + qA (x, t) (3.6) 

c i = i 

and the coefficients in the equations (3.4) are 

F iai ... ak (x) = 0, fc 7^ 1, i = 1,2,3, ai, ...,a fc = 0, ...,3, (3.7) 



Fa(lH ^_I^), i=1 , 2 ,3. 
or 1 c at 



(3.8) 



We define 
Then the identity 



F 00 = 0, F 0i = -F i0 , i = 1,2,3. (3.9) 
3 dtdx a dt dx? 

f (3 - 10) 

similar to the identity (3.3) holds. By making use of the second identity (3.2) and the 
relations (3.8) - (3.10) we can rewrite the equation (3.4) with the coefficients (3.7), (3.8) as 
the equation with Lorentz force 

dt d (dtdx»\ q ^ 3 dtdx v 

mC TsJt \ds~dt) ~ ~~c V ho v ds~dt' 

, , <9A,(x, t) dAJxA) 

= _ „,„ = 0,...,3. (3.11) 

The interaction is defined by the product of the charge g and the vector potential A M (x, £) 
of the external field. 
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The relativistic Coulomb law is given by Lorentz covariant equations describing the elec- 
tromagnetic interaction of two particles with the charges qk, k = 1,2, 



dt d ( dt dx£\ q k ^ 



x 



c- 1 (l-c- 2 


dx k 




dt 



dt dx v k 
j) ds k dt ' 

-1/2 



Fj;Hu(Xk, Xj) — 



dt 
ds k 

dAj-v^Xk, Xj) dAj-^Xk, Xj) 



dx\ 



(3.12) 
(3.13) 



for any permutation k,j of the integers 1,2. The world lines Xk(t), k = 1,2, satisfy the 
condition x° k (t) = ct where c is the speed of light. Lienard - Wiechert vector potentials are 
given by the following relations 



3 

Aj.^Xk, Xj) = -AnqjK ^ rj^ / dte {x k - Xj(t)) 

, — n J 



q j Kr] i 



d 



u = 
3 



dx v 3 (t) 
dt 



d 



-Xa 



f//(() pW))J ^|x fe -x^t(0))|-^(4-4(t(0)))^4W0))J , 

x2-ct(0) = |x fc -x J -(*(0))| (3.14) 

where e (x) is the first distribution (2.44). By making change of the integration variable 

xl-ct(r) = (\x k - X Mr))\ 2 + r) l/2 (3.15) 

it is easy to prove the relation (3.14). For r = the relation (3.15) coincides with the last 
relation (3.14). 

For a world line Xj(t) we define the vector 



dt5(x — Xj(t)) 



dXj(t) 
dt 



(^(c-V))5(x-x,(c-V)), M = 0,...,3. 



(3.16) 



The world line Xj (t) satisfies the condition x®{t) = ct. Hence the definition (3.16) implies 
the continuity equation 



d 

dx° 



A*, Xi) = - ± ( f y)) JLs ( x - X) („->*»)) , 



d T i( \ _ ( u % 

— J [X,Xj) - \ ~^ x - 



dx 



4 'c^x 



-£-5(x-Xj(c- l x°)),i = l,2,3, 



(3.17) 



The world line Xj(t), x° = ct, is called time like, if it satisfies the condition 



rfx j(*) 
dt 



< c, 
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For the time like world line xf(t) the supports of the distributions (3.16) lie in the closed 
upper light cone and 

/ dte (x - = / Aeo(^ -y)J dtS(y - x 3 (3-18) 

The supports of both sides of the equality (3.18) lie in the closed upper light cone. The 
distribution e (x) satisfies the equation (2.42). Hence both sides of the equality (3.18) 
satisfy the equation 

(-(d x ,d x ))nx) = J^x,x J ). 

The difference of two solutions of this equation is a solution of the equation (2.42). Since the 
solution of this equation similar to the solution of the equation (2.42) is unique in the class 
of the distributions with supports in the closed upper light cone both sides of the equality 
(3.18) coincide. 

The relations (3.17), (3.18) imply the gauge condition for the vector potential (3.14) 

T,rr-^A j . ll (x 1 x j ) = 0. (3.19) 

/i = 

Due to relation (3.19) the tensor of the electromagnetic field (3.13), (3.14) satisfies Maxwell 
equations with the current 4irqjK7] fJil j i J fl (x k ,Xj). 

In contrast to the equation (1.3) the equations (3.12) - (3.14) satisfy the causality con- 
dition. The vector potential (3.14) depends only on the world line points Xj(t) lying in the 
closed lower light cone with the vertex at the point x k . 

In the equations (3.12) - (3.14) the electromagnetic interaction is defined by Lienard - 
Wiechert vector potentials. Let us consider the interaction coefficients 

A$\x) = K jk . Ql v (d x )e ofi (x) + K^Q^dJeoflix), 

QlM) = ~V,u(d x , d x ), Ql v {d x ) = A q^q-^- Vr@x, d x ). (3.20) 

The interaction coefficient (3.20) is defined by two constants Kj k . , Kj k -i. Due to the re- 
lation (2.45) for Kjk ; i = the interaction coefficient (3.20) is equal to r] fll/ Kj k . : oe (x). The 
polynomials Q l ^ u (d x ), I — 0, 1, are connected to Minkowski metric 

V,u E V a(3 Qlp(d x )=4QUd x ), E V aP QUd x )=0. (3.21) 

o,/3 = 0,...,3 o,/3 = 0,...,3 

We change in the equations (3.12) - (3.14) the interaction coefficients A.nq k q,jC~ x Ki]^ v e (x) 
for the interaction coefficients (3.20). The relation (3.18) is still valid when the distribution 
eo,o(x) substitutes for the distribution eo(x). By making use of the obtained relation and 
the relations (2.45), (3.17), (3.20) it is easy to prove the relation 



3 r dx u (t) 3 r dx v (t) 

E / dtAf\x k - xj(t))-^ = (K jk , + K jkil ) E V,u / dte (x k - xj(t))-^. (3.22) 



The constants Kj k -i, I = 0, 1, of the interaction coefficient (3.20) are included in the vector 
potential (3.22) as the sum Kj k . Q + Kj k .\. 
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Let us introduce the interaction coefficients into the equation (2.41). Let j, k be the 
permutation of the numbers 1,2. We construct the equation for j particle. Let us multiply 
the equations (2.41) for the particles 1 and 2. We change the differential operator 



-^H-^) (3 23) 



for the differential operator 

( . d \ ( . d 



dx v i I \ dx 



+ Ai%\(x J -x k ). (3.24) 



(We do not change in the equation (2.41) the differential operator —id v for the differential 
operator — id v + qA u (x)). In the differential operators (3.24) the interaction coefficients 
A^ k J (x) are given by the relations (3.20). For any matrix A E SL(2, C) the equality 



'3 \ 3 



€ EAj(ry = £ Kl(A)A%(A)AW 2 (x) (3.25) 

\A t = / Mi)M2=0 

holds. The differential operator (3.24) transformation is similar to the differential operator 
(3.23) transformation. The support of the distribution (3.20) lies in the closed upper light 
cone. Hence the differential operator (3.24) differs from the differential operator (3.23) only 
for the vector Xj — x k lying in the closed upper light cone. In order to obtain an equation of 
the type (1.2) we integrate our equation with respect to the variable Xk 

E E / d 4 x k 

n s = l,...,4Z s +2, s = 1,2 ii s = l,...,2/ s +l, s = 1,2 

m s rh s ,n a h a ] (fs)n s n s (%s) ~l~ 

2 

II (7?. ) f.(^))m.m. ) n.n.*(a;2)(/.)„ i n.(+0,X 8 ) + 



vs = l \u = 



s = l 



3 



^ - ^) II (^.(^))m.*.,n.ft.(/.)n.ft.(x.) = 0. (3.26) 

1^1,^2 = S = l / 

In the equations (3.14) we integrate along the particle world line the product of the interac- 
tion coefficient and the velocity vector 



3 

-AirqjK ^ / dte (x k - Xj(t)) 



dx»{t) 
dt 



(Due to the relation (3.22) the interaction coefficient (3.20) gives a similar interaction term.) 
A particle has no a world line in the quantum mechanics. The particle probability density is 
given by the wave function. In the equations (3.26) we integrate over the space the product 
of the same interaction coefficient and the vector constructed by means of the wave function 



/ 



3 4/ fc +2 2/ fe +l 

d Xk A^J\xj — Xk) ^ ^ (ll k j k (A))m k m k ,n k h k {fk)n k n k {Xk)- 

u = n k = \h k = \ 
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The first and the second terms in the left - hand side of the equation (3.26) correspond to 
the multiplied equations (2.41) for the particles 1 and 2. The equation (3.26) transforms like 
two equations (2.41). 

The relations (3.20) define the interaction coefficients A^J) 2 (x) in the equations (3.26). 
The degree of the homogeneous polynomials Qi ll/2 (d x ), I — 0, 1, is 21. Hence the second 
relation (2.44) implies that the supports of the distributions Q l (d x )eo,o(x), I = 0, 1, lie on 
upper light cone boundary. Thus the equations (3.20), (3.26) satisfy the causality condition. 
The interaction propagates at the speed of light. 

The definitions (2.49), (3.20) imply 

A^\x) = \lm(K krfi + K krA )(2ny 4 J d 4 pexp{-i(p, x)} x 

« p0 + ^ - CX^-iPij ■ (3 - 27) 

Due to ([3], Chapter 12, relation (12.91); [12], Chapter III, relation (3.58)) the propagation 
function of the vector particles in the Yang - Mills theory is 

Df v {x) = hm-5 ab (27r)- 4 J d 4 pexp{-*M} x 

((p,p) + ie)' 1 - (1 - ") fop^ e ) • ( 3 -28) 

The number a is the consequence of the gauge condition for the value d^A 11 . The choices 
a — 1 and a = are called the gauge conditions of Feynman and Landau ([3], Section 12.2.2). 
The distribution (3.28) differs from the distribution (3.27) in the rule of going around the 
poles in the integral. 

We have considered up to now that the interaction propagates at the speed of light. 
Let the interaction propagate at the speed less or equal to the speed of light. Hence the 
interaction coefficients (3.20) may be changed in the following way 

A^\x) = ]T / dX 1 dX 2 K kr ^X 1 ,X 2 )Ql u (d x )exlxl(x). (3.29) 
i = o J 

If the distributions K k j.i(X 1: A 2 ) = K kj .i5(X 1 )5(X 2 ), I = 1,2, then the interaction coefficient 
(3.29) coincides with the interaction coefficient (3.20). The interaction coefficients (3.29) 
satisfy the covariance relation (3.25). Changing in the equations (3.26) the interaction 
coefficients (3.20) for the interaction coefficients (3.29) we obtain in the general case the 
theory where the particles interact between each other by means of the massive particles. 

We consider the equations (3.26), (3.20) for the case h — l 2 — 1/2, l x — l 2 — 0. We 
suppose that the constants K 12 -i, 1 — 0,1, are small. If we neglect the term containing 
the interaction coefficient A^J 2 (x), then the functions {f2)m 2 { x 2) are the solutions of the 
equation (2.41). It is possible simply to assume K\ 2 -i — 0, I — 0, 1. We do not know the 
physical reason of the assumption Ki 2 -i = 0, K 2 \-i ^ 0, I — 0, 1. However the differential 
equations (3.20), (3.26) with these constants are mathematically self - consistent. 

Let the functions (f 2 ) m2 (x 2 ) be the solutions of the equation (2.41) and be given by the 
relation (2.51). Let the integral 

^(/ 2 v( + o,y 2) = { J /^ ( ,fritn y2) ; '" 2 = \ 2 ; < 3 - 3o > 

I / d 3 y 2 (f 2 ) m2 - 2 (+0, y 2 ), m 2 = 3, 4 



/ 
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be not equal to zero for some index m 2 . The equations (2.41), (3.26) imply 

E E(t") 

mini I + ^ 21) (^l) (fl)m(xi) + 

ni = lu = \ OX l / 

£ ()9(^?))m 1 n 1 (/i)n 1 (a:i) = (A°U (+0, Xl ) , (3.31) 

m = 1 

£( 21) (*i) = i (/ rf 3 y 2 (/ 2 ) m2 (+0,y 2 )) _i x 
]T / ^^(X! -x 2 ) £ (T^Un^kO^). (3.32) 

fl = U2 = 1 

The equation (3.31) is similar to the equation (1.2). The right - hand side of the equation 
(3.31) defines Cauchy problem. The first particle interacts with the field of second particle. 
The vector potential of this field is given by the equality (3.32). The substitution of the first 
relation (2.44) and the relation (3.20), K 2 i-i = 0, into the vector potential (3.32) yields the 
delayed vector potential of second particle electromagnetic field 



)m2ri2 ^ 



B^\x 1 )=(fd"y 2 (f 2 ) m2 (+0,y 2 )Y £ £1^(7* 

VJ ' i/i = 0m2 = 1 

(Ani)- 1 K 21 . J d 3 x 2 |x! - x 2 | _1 (/ 2 )n2(^? - |*i -x 2 |,x 2 ). (3.33) 

If we choose in the vector potential (3.33) the functions (f 2 ) m2 (x 2 ) = (/2)m 2 ^( x 2), we get 
the vector potential of Coulomb type. The second particle moves freely and the functions 
{h)m 2 { x 2) are the solutions of the equation (2.41) given by the relation (2.51). It is possible 
to choose the initial functions (/2)m 2 (+0, x 2 ) only. 
The following equality 

/ d 4 ye^ ^ n i {x - y)e^(y) = e^^(x). (3.34) 

holds. Both sides of the equality (3.34) have supports in the closed upper light cone and 
satisfy the equation (2.42). The solution of the equation (2.42) is unique in the class of the 
distributions with supports in the closed upper light cone. By making use of the relations 
(2.51), (3.20), (3.34) it is possible to rewrite the equality (3.32) as 



B^\x) = (/ rf 3 y 2 (/ 2 °) m 2(+0,y 2 )) _1 



x 



E E ( E (l Ul l U2 )m 2 n 2 (-i^) - (7^2))™ 



X 



vi = ri2 = 1 \ v 2 = 



/ d A y{K 21]0 Ql Vi {d x ) + K 21 ., 1 Ql^(d x ))e 0fiA (x-y)S(y ){f^) na (+0,y). (3.35) 

The relations (2.47), (2.48), (3.34) imply the equality 

00 

e* fiA {x) = 0{7?)0{{x,x)) J2^ 2 n (-ir +1 (2n4 n+2 (n + l)\(n + 2)\r\x,xr + \ (3.36) 



n = 
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By making use of the relations (2.21), (3.20), (3.36) we can rewrite in the open upper light 
cone the vector potential (3.35) with the functions (/2)n 2 (+0, y) = (/2)n 2 ^(y) as 

B^\x) = + B^ l \x) (3.37) 

where the constant vector potential 

£f i0) = -(Stt)" 1 ^ ( E ( 7 °Un 2 (/ 2 ) J E ^( 7 V/9(^)) m2 „ 2 (/ 2 )„ 2 (3.38) 

\n 2 = 1 / n 2 = 1 

and the remainder vector potential 

B? 1;1) (z) = 

/ rf 3 y 2 (/ 2 °) m2 (+0,y 2 )) ^r(-ir 1 (2rf +2 (n + l)!(n + 2)!)- 1 x 



E E (E ^V 2 /)™ ( 

v\ = n 2 = 1 \f 2 =0 \ 



(/ 2 )n 2 (^21;0Q°, 1 (^)+X 21;1 QL 1 (^))(x,x)" +1 . (3.39) 

It is possible to neglect the vector potential (3.39) in the right - hand side of the equality 
(3.37) for the small /i 2 . For /i 2 — the vector potential (3.37) is simply equal to the vector 
potential (3.38). 

The equation (3.31) is written for the functions (fi) ni (xi) equal to zero for the negative 
x\. Let us consider in the open upper light cone the equation of the type (3.31) for the 
arbitrary functions (tpi) mi (xi) and the vector potential (3.38) 

E E^)™ f-i^ +B^°A (VlWxi) + E (^?))mi»i(^l)n 1 (xi) = 0. (3.40) 
ni = lv = \ OX l J m = l 

The equation (3.40) is similar to the equation (1.2). Let the matrix a{ji\) act on the equation 
(3.40). In view of the relations (2.21) we have 

E EKWm f"<A + Bl 21 'A (6U*i) = 0, (3.41) 

ni = lu = \ OX l ) 



^iU(zi) = E (/5(^r 1 ))mm 1 (6)n 1 (a:i). (3.42) 
m = 1 

Then the solution of the equation (3.40) in the open upper light cone is 

4 oo 

(V>lU(zi)= E E(^0" 1 (-^l) fc (/3(^r 1 )C fc )m 1 n 1 (ei)n 1 (3.43) 
ni =1 fe = 

where a vector(£i) ni is independent of the variable X\ and the matrix 

C m = E ^ 21; ° ) (7°7 l/ )m 1 n 1 + ^(7°)™ (3.44) 

;/ = 
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is Hermitian if all components of the vector potential (3.38) are real. 

The vector potential (3.38) and the mass fi\ define the energy spectrum of the solution 
(3.43). Choosing the vectors (/2)n 2 it is possible to obtain a practically arbitrary energy 
spectrum. 

The vector potential (3.35), (/2)n 2 (+0> x 2) = (/2)n 2 ^( x 2), is equal to zero out of the closed 
upper light cone. The functions (fi) ni ( x i) satisfy the free equation (2.41) out of the closed 
upper light cone. If the initial wave function (/i)m(+0, Xi) is equal to (/i) ni 5(xi) (at the 
initial moment both particles are at the coordinates origin), then in view of the relation (2.51) 
the function (/i)n 1 ni(^i) is equal to zero out of the closed upper light cone. We study the 
equation (3.31) with the singular vector potential (3.35), (/2)n 2 (+0> x 2) — (/2)n 2 ( K x 2) an d 
the singular initial wave function (/i) ni (+0, Xi) = (/ 1 ) ni 5(x 1 ). The solution of the equation 
(3.31) is also singular. The solution singularities lie on the upper light cone boundary. 
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